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The recent global analyses of the nuclear parton distribution functions (nPDFs) lend support to the 
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1. Introduction 

The global analyses of nPDFs approach the experimentally observed differences between the 
hard-process cross-sections involving bound and free nucleons, (j boundnucleon ^ (j free nucleon ^ w j m 
the toolkit of perturbative QCD (pQCD). Different processes in different kinematical regions are 
studied in order to extract the universal — process independent — nuclear PDFs, and to estimate 
the conditions in which the factorization theorem of QCD appears applicable. In this talk we 
recapitulate the essentials of the latest offspring of the EPS -generation |jl[ |[ EPS09 and 
release a piece of news about the recent EPS -activities. 

2. Analysis 

2.1 Experimental Input 

The main source of experimental input in EPS 09 comes from the measurements of deep in- 
elastic scattering (DIS) and Drell-Yan (DY) dilepton production. These data alone, however, leave 
the gluon PDFs badly unconstrained. As noted in [Q], part of such freedom can be reduced — 
in a way that do not cause disagreement with other data — by adding data for average high-/?r 
pion production in d+Au collisions measured at RHIC. It should be noted that the fragmentation 
functions which in pQCD-calculations turn the partons to pions might also experience a modifica- 
tion with respect to the free proton fragmentation functions [§J. This effect is not accounted for in 
EPS09 where the pion-observables are mainly needed for constraining the nuclear gluons. 

2.2 Theoretical Framework 

In EPS09, the bound proton PDFs for each flavor i are linked to the free proton ones by 

y; boundproton (^,G 2 ) =/?f(x,G 2 )/; freeproton (x,e 2 ). (2.1) 

The nuclear modification factors Rf are parametrized at charm quark mass treshold Q 2 = Q 2 = 
1 .69GeV 2 , assuming a flavor-independent nuclear modification for light sea quarks R^ = Rj = Rj, 
and valence quarks Rf ly = Rj at Qq. This is, of course, just a simplifying assumption and our plan 
is to to release this assumption in a future. 

The EPS09 employ the MS, zero-mass variable flavor number scheme. Accordingly, the base- 
line free proton PDF set used was the CTEQ6.1M [|7|], which probably remains as the last CTEQ- 
set in this simple scheme, the more recent ones involving a more complete prescription to treat the 
heavy quarks. Recently, we have also implemented such a scheme to our fitting engine. However, 
the available nuclear data are not very sensitive to the heavy quarks, and our tentative results do not 
display large deviations from the zero-mass ones. Therefore, the usage of EPS09 with newer sets 
of the free proton PDFs like CT10 [JsJ] should be justified. 

2.3 Fit Quality 

In order to find the Rfs that optimally reproduce the experimental data, a certain quality crite- 
rion is need. We minimize a generalized # 2 -function which is essentially 

N ieN \ u > / 
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where D, and a, are the experimental data points with their errors, and 7} is the corresponding 
calculated value. A small extension to this due to the additional normalization uncertainty of the 
pion data is explained in The additional weights wn are assigned to emphasize certain data 
sets with only a small amount of data points but carrying physically relevant information. Partly, 
this is a technical issue aiming to improve the convergence of the fit with more open parameters, 
but certainly introduce some subjectivity to the central fit. Unfortunately, the weights also distort 
the error analysis explained shortly, and inevitably lead to a somewhat underestimated PDF-errors. 
Reducing the impact of such reweighting is one of the issues we aim to improve in the future. 
2.4 Error Analysis 

The uncertainty analysis propagates the experimental errors from the data to the PDFs. The 
Hessian method, used in EPS09, relies on a quadratic approximation of x 2 around its minimum Xq 

X 2 {Wi}) « Zo +ZSa i H ij Sa j = * 2 + £z 2 , ^ 

'j i 

where the 5a,s are the deviations of the original fit parameters a,- from the best-fit values. The 
key point is that the Hessian matrix Hij can be diagonalized thereby finding an uncorrected basis 
{zi}. Relying to so called 90% confidence criterion, we find that that distorting the best fit in the 
directions of each Zi such that the x 2 changes by 50, still gives a resonable agreement with the 
data. These corners of parameters define the uncertainty sets Sf facilitiating the estimation of the 
upper \lower error in any PDF-dependent quantity X with best-fit prediction X(Sq), by 



AX + \- = /£max\min [X(S+) -X(S ),X{Sr) -X(S ),0] 2 . (2.4) 



2.5 Results 

We summarize the resulting nuclear effects in PDFs by Figure [I] which shows the nuclear 
effects in Lead at a low scale Q 2 = 1.69GeV 2 and also at a higher scale Q 2 = lOOGeV 2 to demon- 
strate the effect of DGLAP evolution. For lighter nuclei the effects are closer to unity. The 
blue bands denote the EPS09-uncertainty. We also show the results from two alternative analy- 
ses [12, 13] helping to appreciate the well and badly constrained components of the nuclear effects. 

3. Recent Activities 

For studying the universality of nPDFs further, data from other processes are always welcome. 
One of such processes is the neutrino beam induced DIS. It has been suggested [Q] that the NuTeV 
vFe-data [ |iot ] favors different nuclear effects than extracted from charged lepton DIS and DY 
processes. However, our analysis [[II]] with a larger set of data implies that this is a sheer artefact 
brought about looking only at the NuTeV data that appear to contain internal incosistencies. 

In order to cleanly probe the nPDFs, an electron-ion collider like the planned LHeC would 
be an ideal machine. By fitting to a simulated LHeC data, we have found that such data would 
dramatically improve the determination of small-* gluons and sea quarks. Another promising probe 
for nPDFs, that could be available before the realization of an electron-ion collider, would be the 
proton-ion runs at LHC. We have found that e.g. production of heavy electroweak bosons in such 
a collisions would provide a surprisingly neat tool to investigate the nuclear modifications to the 
PDFs with not much sensitivity to the underlying set of free proton PDFs. Both results mentioned 
here will be published shortly. 
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Figure 1: The nuclear effects for valence quarks Ry, sea quarks R$, and gluons Rq in Lead. The NLO results 
from EPS09, HKN07 iQ, and nDS Q are shown at two values of Q 2 . 
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